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Abstract 

The subject of the charged-to-neutral yield ratio for BB and DD pairs near their 
respective thresholds in e + e~ annihilation is revisited. As previously argued for the B 
mesons, this ratio should exhibit a substantial variation across the T(45") resonance due 
to interference of the resonance scattering phase with the Coulomb interaction between 
the charged mesons. A simple alternative derivation of the expression describing this 
effect is presented here, and the analysis is extended to include the D meson production 
in the region of the ip^nty resonance. The available data on Kaon production at the 
0(1020) resonance are also discussed in connection with the expected variation of the 
charged-to-neutral yield ratio. 



1 Introduction 



The near-threshold resonances 0(1020), 0(3770), and T(4S), decaying respectively to pairs 
of pseudoscalar mesons KK, DD, and BB are well known 'factories' for production of these 
mesons in e + e~ annihilation. In a number of experimental studies it is important to know 
the relative yield of the pairs of charged and neutral mesons in the corresponding resonance 
region, i.e. the ratio 

Rc/n _ *(e+e- -» P+P-) 

a(e+e- -> P°P°) ' { ) 

where P stands for the pseudoscalar meson, i.e. K, D, or B. The knowledge of this ratio is 
of particular importance for the studies of the B mesons at the P-factories, and dedicated 
measurements have been done at the T(4S) resonance by CLEOpjJE], BABAR[3 t 0] and 
Belle jH]. The central values of the ratio R c ^ n found in these measurements typically range 
from 1.01 to 1.10 with the latest result @] being 1.006 ± 0.036(stat) ± 0.031 (ays). The same 
ratio for the production of D meson pairs is likely to be studied in detail at the -0(3770) 
resonance by the imminent CLEO-c experiment [6 . For the Kaon pair production the most 
detailed available scan of production of charged and neutral mesons in the resonance region 
has been done by the SND collaboration 

The theoretical treatment of the ratio P C//n at these three thresholds near the corre- 
sponding resonances bears considerable similarity, although specific factors contributing to 
the ratio in each case are slightly different. In a way, the simplest case is presented by the 
threshold B pair production, where the deviation of P C//n from one is essentially entirely 
determined by the Coulomb interaction between the charged B mesons, since the mass dif- 
ference between the B mesons is very small: m^o — m B + = 0.33 ± 0.28MeV[B]. For the 
DP) pair production at the 0(3770) peak the difference in the P-wave kinematical factor 
p 3 due to the substantial mass difference between the charged and neutral D mesons is the 
main source of deviation of R c l n from one, and the Coulomb interaction effect is somewhat 
smaller, but is still of a measurable O(10%) magnitude. Finally, in the 0(1020) resonance 
region in addition to the mass difference and the Coulomb interaction effects, the amplitudes 
of production of K + K~ and K^Kg also differ due to a non-resonant isovector, 1=1, con- 
tribution coming from the electromagnetic current of the u and d quarks. Another technical 
difference between the heavy mesons and the Kaons is that the former can safely be consid- 
ered nonrelativistic at energies in the region of the corresponding resonance (the velocities of 
the mesons at the corresponding resonances are v(B)/c ~ 0.06 and v{D + )/c ~ 0.13), while 
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the velocity of the charged Kaons at the <fi resonance is v(K + )/c ~ 0.25, and the 0(v 2 /c 2 ) 
relativistic effects can show up at some level of intended accuracy in R c l n . 

The kinematical effect of the mass difference, and of the non-resonant isovector amplitude 
for the case of Kaon production, lead to factors in the R c ^ n , which are rather smoothly varying 
within the width of the corresponding resonance. The behavior of the Coulomb interaction 
effect is however quite different. Namely, as recently pointed out [5], with a proper treatment 
of the strong resonant interaction between the mesons, there is a phase interference between 
the Coulomb interaction and the resonance (Breit-Wigner) scattering phase. Since the latter 
phase changes by ix across the resonance, the effect of the Coulomb interaction in R c l n 
should exhibit a substantial variation with energy within the width of the resonance. Thus 
in particular a comparison between different measurements of the charged-to-neutral ratio at 
the T(45) peak can only be meaningful with proper account for the differences in the specific 
experimental setups such as the beam energy spread, stability, calibration, etc. Needless to 
mention that the best approach would be a dedicated scan of the T(4S) peak. The details of 
this variation depend jH| on the specifics of the transition from the strong interaction region 
at short distances to the long-range Coulomb interaction and also on such details as the non- 
resonant scattering phase for the mesons. A dedicated scan of the charged-to-neutral ratio 
thus could possibly provide an information on rather fine properties of strong interaction 
between the heavy mesons, which are not likely to be accessible by other means. It may 
well be that a more experimentally feasible object for such study is provided by the ^(3770) 
resonance, and possibly can be explored in detail in the forthcoming CLEO-c experiment. 
The main motivation of the present paper is to extend the analysis of Ref . jU] to the case of 
the DD production in the region of this resonance, and to estimate the magnitude of the 
expected effects, in the extent that it can be done within the present (rather approximate 
at best) understanding of the details of the strong interaction between heavy mesons. 

The typical amplitude of the varying part of the ratio R c ^ n due to the Coulomb interaction 
is set by the Coulomb parameter ct/v for the charged mesons, which tells us that for the 
D mesons at ^(3770) the amplitude of the variation should be approximately twice smaller 
than for the B mesons at T(4S) (and another factor of two smaller for the Kaons at 0). In 
either case the Coulomb parameter is sufficiently small to justify limiting the consideration to 
only the linear term in the Coulomb interaction treated as perturbation. A simple derivation 
of the linear in the Coulomb interaction correction to production of a charged meson pair 
with the exact treatment of the strong interaction scattering is presented in Section 2. This 
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derivation is slightly different from that given in Ref.jH], and hopefully, more transparently 
leads to the expression for the 0(a/v) term found there. In Section 3 are presented specific 
estimates of the behavior of the ratio R c / n across the -0(3770) resonance, which take into 
account both the mass difference between the D mesons and the effect of the Coulomb 
interaction. 

It can be noticed that the very fact of the existence of the variation of R c l n within the 
resonance width is essentially model independent and is a direct consequence of a somewhat 
standard application of the quantum mechanical scattering theory. Nevertheless, there has 
been some scepticism expressed regarding the existence of the predicted jU] rapid variation 
of R c/n at the resonance, most notably in the introductory part of the recent experimental 
paper [4,.. Although the reported in Ref. [4J measurement has been done at just one point in 
the energy at the peak of the T(45) resonance, the paper implicitly puts the prediction of the 
variation in doubt by arguing that "such rapid variation in the charged-to-neutral ratio has 
not been observed in scans across the 0(1020) resonance" (with a reference to the scans by 
the SND collaboration^). Thus it appears to be useful to point out the difference between 
the predictions for the B production at T(4S) and the Kaon production at 0(1020) following 
from the discussed approach. The main difference is that the amplitude of the variation at 
the is expected to be suppressed by a factor of approximately 0.25 due to about four 
times larger velocity of the Kaons at the relevant energy, and additional differences in the 
details arise through different parameters (e.g. widths) of the and T(4S) resonances. A 
comparison of a 'representative' theoretical curve with the behavior of R c l n across the 
resonance extracted from the scan data[7j is presented in Section 4. Given the existing 
theoretical and experimental uncertainties, a meaningful comparison can only be done at a 
mostly qualitative level, and it is left entirely up to the reader to assess from the presented 
comparison whether the scan data indeed exclude a variation with expected amplitude of 
the charged-to-neutral ratio across the resonance. 



2 Coulomb interaction correction to production of charged 
mesons 

If the mesons were point particles devoid of strong interaction and produced by a point 
source, the Coulomb interaction between charged mesons would enhance their production 
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cross section by the textbook factor 



F . = i + ™ +0 (<»v) , 



(2) 



where v is the velocity of each of the mesons in the center of mass frame, which was the early 
prediction PU] for the ratio R c ^ n at the T(AS) resonance. It was however then realized that 
the effects of the meson electromagnetic form factor JI] and of the form factor in the decay 
vertex T(4S) — > -B-B[l2J|H] tend to reduce the enhancement. In all these approaches the 
Coulomb rescattering of the produced mesons is considered, assuming that the propagation 
of the meson pair between the production vertex and the rescattering as well as after the 
Coulomb rescattering is described by a free propagator, and thus ignoring the fact that the 
strong interaction introduces the scattering phase factors exp(i5) in these propagators. This 
would be a reasonable approximation if the strong phase 5 were small. In the region of a 
strong resonance however this assumption is definitely invalid, since the phase changes by 7r 
within the width of the resonance. The calculation described in Ref.[9 allows to completely 
take into account the strong scattering phase as well as the Coulomb interaction between 
the mesons. In the first order in the Coulomb rescattering the expression [Hj for the Coulomb 
correction factor has a simple form 



where p is the meson momentum and V(r) is the potential for the rescattering interaction. 
Clearly, up to a possible form factor, V(r) = —a/r for the Coulomb interaction. The 
short- distance cutoff parameter a in eq.Q accounts for the fact that in the region of a 
strong interaction at short distances the mesons spatially overlap and in fact the system 
is a mixture of heavy and light quarks and gluons rather than a state of two individual 
mesons. Therefore at such distances there is no separation between the states with charged 
and neutral mesons and thus any difference in their Coulomb interaction disappears. If one 
introduces form factors in the Coulomb interaction in order to account for the structure of 
the mesons, the interaction potential V(r) would automatically vanish at short distances 
and the ultraviolet cutoff in the integral in eq. (JHJ) would be provided by the form factors. If 
as in the earlier analyses, one sets the strong scattering phase 5 to zero, the expression in 
eq.(JH} is finite even if an unmodified Coulomb potential is taken down to a = 0, in which 
limit it reproduces eq.(J2J). However at any finite S there is an essential dependence on the 
ultraviolet cutoff parameter a which cannot be removed. 
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In this section we present a slightly different from that in Ref. |5J derivation of eq.® 
by considering the following problem. Two mesons, each with mass m interact strongly at 
distances r < a. The strong interaction in the P wave at the total kinetic energy E produces 
the scattering phase 5. At longer distances, r > a the mesons interact via the potential 
V(r). The meson pair is produced by a source localized at distances shorter than a (this 
obviously corresponds to the production in the e + e~ annihilation). The problem is to find 
the correction of the first order in V to the production rate. 

In order to solve the formulated problem we consider the radial part of the wave function 
of a stationary P-wave scattering state, written in the form R(r) = x{ r )/ r - According to 
the general scattering theory [Ti] the asymptotic form of x{ r ) at r ^ oo in the absence of 
the long-range potential V(r) is 

X {r) = 2 cosGor + 5) = e lS e ipr + e^ 5 e~ ipr . (4) 

Furthermore, in the absence of the potential V the function x(r) satisfies the Schrodinger 
equation for free motion in the P wave: 

X "(r) + (V - 1) X (r) = (5) 

at all distances r outside the region of strong interaction, i.e. at r > a. Thus the solution 
valid at all such distances, and having the asymptotic form (j3J) can be written explicitly: 

X(r) = e l5 f(pr) + e- l5 r(pr) , (6) 

in terms of the function 

and of its complex conjugate /*. At distances r < a the dynamics is unknown to the extent 
that continuing description of the system in terms of the meson pair wave function in the 
region of strong interaction does not make much sense, since most likely such description 
should involve entirely different degrees of freedom. Nevertheless whatever complicated the 
'inner' dynamics of the system may be, the wave function in eq.((5|) at r = a provides the 
boundary condition for the 'inner' problem. In particular, it determines the normalization of 
the state for the 'inner' problem. Thus the rate of production of the system, and consequently 
of the meson pair, by a source localized at r a is proportional to |x(a)| 2 . 

When the long-range potential V(r) is turned on and considered as a perturbation the 
wave function at r > a changes to x( r ) + $x{ r )- Once the corrected solution is normalized 
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at r — > oo in the standard way, i.e. as the asymptotic wave in eq.(JH), the normalization at 
r = a generally changes and hence also changes the production rate. In the first order in V 
the correction factor for the production rate is obviously given by 

fo=l + 2^>. (8) 

The perturbation 5\ of the scattering state wave function can be found in the standard 
way by writing it as a sum of outgoing and incoming waves: Sx{r) = 5x+{r) + <5x_(r) , 
similarly to eq. (J5|) . and where Sx+{r) contains at r — > oo only an outgoing wave (exp(zpr)) 
while 5x~ = $X*+ contains at asymptotic r only an incoming wave. The function 5x+{r) is a 
perturbation of the outgoing wave part of the function x(r) in eq.© and is thus determined 
by the equation 

S X 'i(r) + (p 2 - ^) 5 X+ {r)=m V{r)e iS f{pr) . (9) 

The solution to this equation is found by the standard method using the Green's function 
G + (r,r') satisfying the equation 

(£t + p 2 - ^) G +( r > r ') = 5 ^ r - r ') . ( io ) 

and the condition that G + (r, r') contains only an outgoing wave when either of its arguments 
goes to infinity. The Green's function is constructed from two solutions of the homogeneous 
equation, i.e. from the functions f(pr) and f*(pr), as 

G+(r, r') = -L [f(pr) f*(pr') 6{r - r') + f{pr') f*(pr) 6{r' - r)] , (11) 
zip 

where 9 is the standard unit step function. The solution to the equation Q is then found as 

S x+ ( r ) =me t5 / G+(r, r') V(r') f{pr') dr' . (12) 

J a 

It can be noted that at no point in this consideration a knowledge of the dynamics at r < a 
is required. In particular the integral in eq. ()12}) runs from the lower limit at a, since by our 
assumptions the perturbation potential has support only at r > a. 

It is further important that adding the found perturbation 5x to the wave function does 
not change the normalization of the wave function at r — > oo. Indeed in this limit one has 
r > r' in the Green's function in eq. ()12|) . so that the asymptotic behavior of 5x+{r) should 
be derived from the expression 

Sx + (r) I™ = --e l5 f{pr) ^ V(r') \f(pr')\ 2 dr' , (13) 
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which gives the complex phase of 5% + (r) at the asymptotic distances manifestly orthogonal 
to that of the outgoing- wave part of x( r )> since the integral is explicitly real 1 . In other 
words, at r — * oo the correction changes only the scattering phase by adding the 'Coulomb' 
phase to 5. 

At r = a the perturbation Sx+{o) is found from the same expression (|12|). In this case 
one has r < r' in the entire integration region, so that 

S X+ (a) = --e* s r(pa) ^ V(r') [/(pr')] 2 dr' . (14) 

Using this expression one readily finds the change in the normalization of the wave function 
at r = a due to the potential V, and thus arrives at the expression (J3J). 

As is clear from the presented derivation, the result in eq.Q is valid at arbitrary P-wave 
scattering phase 5, and also generally does not assume any condition for the value of the 
product pa. At energy E near a strong resonance dominating the dynamics of the meson 
pair, the scattering phase is determined by the Breit-Wigner formula 

S = 5 BW + 5 l , e 2 ^- = ^— ^, (15) 

A + ry 

where A = E — Eq is the distance in energy to the nominal position Eq of the resonance, 
and 7 is generally a function of energy such that its value at Eq is related to the nominal 
width r of the resonance as 7(Po) = Y/2. Finally, 5\ in eq. ()15|) is the non-resonant P-wave 
scattering phase. 

The properties of the discussed resonances, T(4S") and -?/>(3770) are determined by the 
two meson-pair channels: P + P~ and P°P . The analytical properties of the scattering 
amplitude in the P wave require that the contribution of each channel to 7 and 5± starts at 
the corresponding meson pair threshold as cubic in the corresponding momentum. Thus at 
the energy close to both thresholds one can parametrize these quantities as 

7 (P) = c x (p 3 + + pi) , 5 1 (E) = c 2 (p% + pi) , (16) 

where p + and po are the momenta of respectively the charged and the neutral mesons at the 

energy E. Due to the very small mass difference between the B mesons, the thresholds for 

1 A minor technical point is that formally the integral is logarithmically divergent at r' — > 00 for the 
Coulomb potential, which might put into question the applicability of the condition r > r' in the entire 
integration region. This is the standard infrared divergence of the Coulomb scattering phase, and it can be 
easily dealt with by temporarily introducing a small photon mass A, so that the potential is cut off by the 
factor exp(— Xr) and the integral is convergent. In the final result one can obviously take the limit A — ► 0. 
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the charged and neutral B mesons practically coincide and so do the momenta p + and po- 
However for the DD pairs the D + D~ threshold is by 9.6 MeV higher than that for D°D°, 
which difference is substantial at the energy of the ^(3770) resonance. 

It should be noticed that neglecting the higher terms of expansion in the momenta for the 
parameters 7 and S\ is justified only inasmuch as the momenta are small in the scale of the 
size of the strong interaction region a s . In the calculation of the Coulomb interaction effect 
it is essential that the distance parameter a for the onset of the Coulomb interaction is not 
smaller than a s : a > a s . Otherwise no restriction on the value of the product pa is implied. 
In Ref. jH] these parameters were reasonably assumed to be approximately equal. However 
it should be emphasized that generally the applicability of the formula (jSJ) for the Coulomb 
correction is separate from the applicability of the first terms of the threshold expansion in 
eq.(fTB|). and these two issues can be studied separately in experiments 2 . 

In conclusion of this section we write the explicit formula for the correction described by 
eq.(|3*|) in the resonance region in terms of real quantities [9J: 

"A 2 - T 2 . A . n . nr . 2 7 A 



F c = l + - 

v 



(A cos25i + B sin25i) - — - — - (B cos2o\ - A sin2o"i) 



(17) 



A 2 + 7 2 v ' A 2 + 7 2 

The coefficients A and B are found as the imaginary and the real parts of the integral in 
eq.(jnj) with the Coulomb potential: 



A 
B 

where 



pa 



1 \ 2 cos 2m 



pa 



1 sin 2u + 

u J u 

2 sin 2m / 1 

— ( 1 ) cos 2m 



u \ u 



du it cos2pa sin2pa 

u 2 pa 2(pa) 2 

du cos2»a sin2»a _ , 

H Ci(2pa) , 



u 2 (pa) 2 pa 
dt „ , s r°° dt 



r,-, \ f z . dt , f°° dt 

bi{z) = suit — and Ci\z) = — / cost — 

Jot J z t 



3 Estimates of the expected variation of R c ^ n across the 
^(3770) resonance 

Here we apply the formulas of the previous section, specifically the equations (|17p. ()18jl. and 

(11611 . to an estimate of the scale of the variation in the ratio R c l n for the DI) pair production 

2 In fact an attempt at detecting higher than p 3 terms in the width parameter for the T(4S) resonance 
has been done by ARGUS collaboration|T3j. However the deviation from the cubic behavior turned out to 
be too small within the experimental accuracy. 
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at the f/>(3770) resonance. In doing so it should be clearly understood that at present one can 
only guess (in some 'reasonable' range) the appropriate values of the scattering phase 5i and 
the cutoff parameter a for the Coulomb interaction (or more generally make a guess about 
an appropriate model for the onset at short distances of the Coulomb interaction between 
the charged D mesons). 

As different from the case of the B meson production at the T(AS), in the charged-to- 
neutral ratio for the D meson production at ?/>(3770) the Coulomb effect multiplies the ratio 
of the p 3 factors: 

R c/n = F c (^j . (19) 

Also, according to the Particle Data Tables [S], the width of -0(3770) is somewhat larger than 
that of the T(45): r«) = 23.6 ± 2.7 MeV as opposed to 14 ± 5 MeV, which tends to 
smoothen the variation of R c l n near the ^(3770) peak. 



3770 
E MeV 




10.576 10.580 10.584 

E GeV 



Figure 1: The energy dependence of the ratio R c ^ n for D and B meson pair production 
in the region of the respective resonance: 0(3770) and T(45) (the center positions are 
assumed at E = 3.770 GeV and E = 10.580 GeV respectively ) for some values of a and 
SxiEo): a- 1 = 200 MeV, S^Eq) = 30° (solid), a' 1 = 200 MeV, S^Eq) = -30° (dashed), 
cr 1 = 300 MeV, S^Eq) = 30° (dashdot), and a" 1 = 300 MeV, 6x{E ) = -30° (dotted). 



A sample expected behavior of the ratio R c / n for the D pair production near the ■0(3770) 
resonance is shown in Fig.l for some 'representative' values of the parameters a and 8\{Eq). 
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For comparison the same curves for the B pairs in the vicinity of the T(4S) peak are also 
shown in a separate plot in Fig.l. It should be noted in connection with this comparison 
that generally it is not expected that the phase 8\ is the same in these two cases, although 
the parameter a viewed as characterizing the electric charge structure in a heavy meson 
is likely to be quite similar for D and B mesons, if one considers both c and b quarks 
as asymptotically heavy. With all the present uncertainty in the knowledge of the strong 
interaction parameters, one can see from this comparison that the expected variation of the 
ratio K c l n is quite less prominent for D mesons at the ?/>(3770) peak than for B mesons at 
the T(45). Hopefully, the amplitude of the variation of R c / n at the ?/>(3770) is still sufficient 
for a study in the upcoming CLEO-c experiment. 



So far the most detailed scan data for the production cross section of charged and neutral 
mesons are available only for the Kaons in the vicinity of the 0(1020) resonancejZj, and it 
is quite natural to discuss whether any hint at the discussed variation of R c / n is indicated 
by those data. A consideration of this ratio at the <fi peak however encounters certain 
peculiarities. As mentioned above, the Coulomb interaction effect is relatively smaller for 
production of K mesons at the resonance, and also the relativistic effects can play a certain 
role. The analysis of the ratio R°' n is further complicated by the fact that the production 
amplitude also receives an isovector contribution, which has opposite sign for K + K~ and 
KlKs and thus changes the charged-to-neutral ratio. (In the vector dominance model this 
contribution is considered as the 'tail' of the p resonance.) Neglecting a smooth variation of 
the non-resonant 1 = 1 amplitude and also of the non-resonant part of the 1 = amplitude 
(e.g. the 'tail' of the uj resonance) one can approximate the formula for R c l n in this case as 



where A\ and A stand for the appropriately normalized relative contribution of the non- 
resonant 1 = 1 and 1 = amplitudes. The charged K mesons are lighter than the neutral 
ones, thus the cube of the ratio of the momenta is larger than one and decreases towards one 
at energies far above the threshold. The datajlj however display a slight general increase 
of the ratio with energy, which shows that the effect of A\ is not negligible. This makes 
any 'absolute' prediction of R c l n in this case rather troublesome, as well as of its general 



4 Discussion of the 0(1020) data and summary 




(20) 
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variation with energy on a scale somewhat larger than the width of <fi. A real fit of the 
involved parameters, including the discussed here rapid variation of the Coulomb factor F c 
requires knowledge of the raw data and of the correlation in the errors. For this reason 
here in Fig. 2 is shown a comparison of the data sets listed in Ref.j7] with a 'representative' 
theoretical curve. 

1 . 6 
1 . 5 



1 . 4 
1 . 3 
1 . 2 

1016 1020 1024 1028 

Figure 2: The SND data[7] at the 0(1020) resonance represented as the charged-to-neutral 
yield ratio for the two reported in JTj scans (PHI9801 and PHI9802) and for two methods 
of identifying the K$ mesons by their decay into neutral or into charged pions: PHI9801 
neutral - stars, PHI9801 charged - diamonds, PHI9802 neutral - triangles, PHI9802 charged 
- squares. The solid curve is explained in the text. 

More specifically, in the paper [7j are listed (in the Table IX) the scan data for the 
cross section for production of separately the K + K~ pairs and the K^Ks pairs acquired in 
two different scans (called PHI9801 and PHI9802). Furthermore the detection of the Kg 
mesons was done by two separate methods, i.e. by detecting their decay into charged pions, 
and into neutral pions. Accordingly the paper [7| lists two separate sets of data for the 
KlK$ production measured by each of these two methods. The systematic error for the 
K + K~ cross section is listed as 7.1%, and for the neutral Kaons as 4% and 4.2% for the 
two identification techniques used. The statistical errors are listed for each individual entry 
for the cross section. The 'data points' in the plot of Fig. 2 are the charged-to-neutral ratia 
calculated from the data listed in Ref . [7] with the errors corresponding to the listed statistical 
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errors only. The 'representative' theoretical curve corresponds to 1.38F C . In other words the 
absolute normalization is set rather arbitrarily, given the described theoretical uncertainty 
in calculating the absolute normalization R c / n at the resonance and given the systematic 
uncertainty in the data. Also any overall variation of the kinematical and the amplitude 
factors in the shown energy interval is neglected. Changing the normalization of the curve 
and also including the overall energy dependence would result in a vertical shift and a slight 
tilt of the curve. The behavior of the Coulomb correction factor in the curve shown in Fig. 2 
is calculated assuming E = 1019.5 MeV, T = 4.26 MeV, cT 1 = 200 MeV, and 8(E ) = 40°. 
The only purpose of the 'theoretical' curve in Fig. 2 is to illustrate the approximate magnitude 
of the expected effect of the variation of R c / n at the 0(1020) resonance, and by no means it 
should be considered detailed prediction. 

Due to large experimental and theoretical uncertainties as explained in the text and can 
be seen from Fig. 2, it is not entirely clear whether the available datajZj exclude or rather 
suggest a variation with the expected magnitude of the ratio R c / n within the width of the 
resonance. Perhaps a detailed global fit of the raw data could provide a statistically 
significant evaluation of the amplitude of such variation. 

In summary. The strong interaction phase significantly modifies the behavior of the 
cross section for production of pairs of charged pseudoscalar mesons near threshold in e + e~ 
annihilation. A simple calculation of this effect in the first order in the Coulomb interaction 
is presented in Section 2. The existence of a strong near-threshold resonance gives rise to a 
rapid variation of the charged-to-neutral yield ratio due to the interference of the resonance 
scattering phase with the Coulomb phase. The specific behavior of this variation is sensitive 
to details of the structure of the mesons and of their strong interaction. Thus measuring 
the discussed effect for B mesons at the T(4S) resonance and for D mesons at the -0(3770) 
resonance might provide an information on these details, which are not accessible by other 
means. The estimated effect for D mesons in the region of the -0(3770) is less prominent 
than for B mesons at the T(4S) but possibly is still measurable. The expected effect of 
the rapid variation of the charged-to-neutral ratio for the Kaon production at the 0(1020) 
resonance is the smallest, and, conservatively, the available data appear to be not conclusive 
enough to confirm or exclude such variation with a statistical significance. 
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